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Abstract We determine the probability distribution of the breaking strength for chains of TV links, 
which have been produced by repeatedly breaking a very long chain. 
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1 Introduction 

Consider a chain assembled from TV random links with independent, identically distributed breaking 
strengths x, which have probability density p(x). It is easy to calculate the probability density p(x\N) 
for the strength of the chain being less than x. In some contexts, however, the relevant question is a 
different one: what are the strengths of chain segments of length TV which are obtained by repeatedly 
breaking a very long chain? Such chain segments are expected to be stronger than their randomly 
assembled counterparts, because they have been produced by a process which has eliminated the 
weakest links. Here we calculate the distribution of the strength of a chain segment of length TV which 
has been produced by repeated breaking a very long chain, of length A/", say. Specifically, we break 
the chain of length A/" at its weakest link, then break each fragment at its weakest link and continue 
the process. We collect the links of length TV produced by this process, and determine the probability 
density of their strengths, p*(x\N). In this paper we obtain formulae which are precise asymptotic 
results for the limit where TV ^> 1. 

The calculation of this distribution involves the use of prior information: we know that the chain 
segment of length TV and breaking strength x was produced by breaking a longer and weaker chain, of 
length TVo and breaking strength xq. What makes the problem difficult is that the prior information, 
in the form of the values of TV and xo, is itself uncertain. It is hard to solve this problem for the case 
of a chain which has previously been broken only once. In the case of a chain which has been broken 
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Fig. 1 Shows p*(x\N) for three different values of N, namely 10 (red), 20 (green) and 40 (blue), for p(x) — 
xexp(-x). Shown are results of computer simulations (histogram) and those of (T| which for the case described 

here takes the form p*(x\N) = N 3 x e~ x {[(l + xf /2 - (1 + x)e x ]e~ 2x + l/2}e~ N( - 1 - (1+x '> e ""> . 

many times, the distributions of xq and Nq themselves depend on previous breakages. The problem 
appears to be very difficult when the chain has previously been broken many times. 

We consider the case where the chain is very long, and we anticipate that there is an asymptotic 
form for the strength distribution which is independent of the initial length A/". We use a self-consistent 
calculation and obtain p*(x\N) in closed form in terms of the cumulative strength distribution of an 
individual link, P(X), and the corresponding probability density p(x) = dP(x)/dx: 



(We use A(N) ~ B(N) to mean that A(N) and B(N) are asymptotically equal in the limit as N — > oo.) 
For comparison, the corresponding probability density for the strength of a randomly assembled chain 
is 



Figure H shows the comparison between equation JQ) and a histogram of data from a numerical simu- 
lation of the chain-breaking process, for a particular choice of link strength distribution p(x). 

We were motivated to study this chain-breaking process because it models fragmentation process 
in which the distribution of strengths of fragments (rather than the distribution of fragment sizes) 
is the principal concern. We are not aware of any earlier investigations of the strength distribution 
of fragments. The chain-breaking process is of interest because it is a fundamental model for such 
processes, which has an explicit solution for the distribution p*(x\N). The model also has direct 
application to the degradation of long chain polymers. 

In section 2 we consider the distribution of the sizes of fragments of a repeatedly broken chain, and 
section 3 describes elementary results for randomly assembled chains. In section 4 these expressions 
are used to obtain a recursion for a sequence of distributions which are related to p*(x\N). Finally 
section 5 shows how the limiting distribution of this sequence is obtained self-consistently, and used to 
derive equation (JJJ above. 

2 Distribution of fragment sizes 

First we discuss the distribution of sizes of chain fragments. After i steps of splitting the chain, we 
have 2 l fragments. Let Wj(iV) be the number of segments of lengths N at step i, and let us consider 
the case where these numbers are so large that it is sufficient to calculate expectation values and ignore 
the statistical fluctuations of Wi(JV) (this assumption is valid in the limit as M — > oo). Using the fact 




(1) 



p{x\N) - Np(x)exp[-NP(x)} . 



(2) 
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that the position of the weakest link has equal probability to be at any site, these numbers satisfy a 
recursion relation 

oo 2 

w l+1 (N)= w=i Wi{M) - (3) 

M=N+1 

Rather than following this iteration for a single chain being broken, it is easier to consider a steady 
state W(N), with destruction of one additional chain being initiated at each step. Thus we seek to 
solve 

oo „ 
M=JV+1 

We find W{N) ~ C/N 2 for N -> oo for some constant C. 



3 Strength of a random chain segment 



By way of preparation we discuss elementary results on the distribution of strengths for a chain with 
completely random links. Let P(x\N) be the probability that a chain of length TV breaks at a tension 
which is less than x. For N = 1, we have P(x\l) = P(x). The probability that a chain of N links 
is unbreakable at tension x is 1 - P(x\N) = [1 - P{x)] N . In the limit where N > 1, the power 
may be approximated by exponentiation: we find P(x\N) <~ 1 — exp[— NP(x)] (which gives J2J) upon 
differentiation) . 

We also require the conditional probability P c (x\xq, N) that a chain of length N breaks at tension 
x if we know that it is definitely not broken by a tension xq. In this case we know that the probability 
of the strength of an individual element being less than x is [P(x) — P(xq)]/[1 — P(xq)], so that 



P c (x\x ,N) = 1 



P(x) - P(x ) 



1 - P(xo) 



N 



provided that x > xq (and zero otherwise). When TV 3> 1, this may be approximated by 

P c (x\x ,N) ~ [1 - exp(-N[P(x) - P(x )])] 0(x - x ) 
where 0{x) is the unit increasing step function (Heaviside function). 



(5) 



(0) 



4 A recursion relation for probabilities of successive breaking tensions 

Our objective is to obtain the probability density p*(x\N) that a chain of length N, which has been 
produced by repeatedly breaking a very long chain, has breaking strength x. 

We determine this distribution of strengths by calculating a related distribution. Consider the 
subdivision of the chain, starting from a very long chain of length TV. When this is repeatedly split, 
for convenience we always discuss the leftmost segment. At the ith stage of subdivision, this segment 
having length N is produced by breaking a link of strength Xi. Let the probability density for the 
strength of this broken link be pi(xi\Ni). The corresponding probability for the ith split to occur at a 
tension less than xi for a segment of length TV is P(xi\Ni). We shall obtain a recursion formula which 
expresses pi+i(xi+i, Ni+i) in terms of this distribution. The chain fragment at stage i + 1 is produced 
by breaking a segment of length iVj, in which all of the links are known to be stronger than a;,-. Both 
of these items of prior information have uncertain values. 

Consider first the distribution of values of Ni, for a given value of 7V i+1 . We have already seen 
(equation that repeated random sub-division of an interval produces a steady-state distribution of 
lengths W(N) ~ 1/N 2 . Subdivision of an interval of length TVj > iVj+i produces an interval of length 
Ni + i with probability 1/iVj. The probability distribution for Ni is therefore proportional to iV~ 3 for 
Ni > N + i (and zero otherwise). For N + i ^> 1, the normalised distribution of N is therefore 

{2N 2 
-Nf N W< N i . (7) 
N +1 > Ni 
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The probability Pi+i(xi+i\Ni+i) is obtained from P c (xi+\ \xi, Ni) by averaging over the probability 
densities of both Ni and Xi. (Note that the latter probability density is related to the unknown function 
that we wish to calculate.): 

P +1 {x l+1 \N l+1 ) = J2 Pn(Ni\N i+1 ) dxi P^Xi+ilx^NMnlNi) . (8) 

N,=N i + 1 + l ^° 

Differentiating JSJ with respect to a^+i and substituting the known expressions for P n (Ni\Ni + i) and 
P c (x\\xo, Nq) we obtain 

00 2N 2 f Xi+1 

p i+1 (x i+1 \N i+1 ) ~ J] — ^± / dx, p(a; i )exp(-iV i [P(a; i+ i) - P(^)])p 4 (^|^) . (9) 

AT=Ar i + 1 + l » "' 

5 Self-consistent solution 

We assume that the distribution pi(x\N) becomes independent of the generation index i when we con- 
sider the subdivision of a very long chain. This leads to a 'self-consistency' condition for the probability 
density pi(x\N) = dPi(x\N)/dx. Approximating the sum in equation by an integral, and replacing 
Pi and Pi+\ by the asymptotic, self-consistent function p^ yields 

Poo (x\N)~2N 2 dN Q -2 / dx exp(-N [P(x)-P(x )])p oo (x \N ) . (10) 
Jn iV o Jo 

We assume that the corresponding probability Poo(x\N) is of the form P oa (x\N) = F(NP(x)) for 
some function F (which increases monotonically from P(0) = to F(oo) = 1). Writing w = P(x), the 
derivative / = F' of F satisfies 

poo \ f W 

f{Nw) =2N dN Q — dw exp[~N (w - w )}f(N w ) . (11) 
Jn iV o Jo 

We express this integral equation in terms of g(X) = exp(X)f(X), and differentiate. We find that the 
primitive of g(X), namely G(X), satisfies 

G"(X) = ^-G'(X) 1g(X) . (12) 

The solution of this equation is of the form 

G(X) = AX 2 + B[exp(-X)(l + X) - X 2 V\{X)] (13) 

where A, B are constants and Ei(Jf) = V J X (y dz e z / z is the exponential integral. The requirement that 
G(0) = implies that B = 0, so that the normalised function f(x) is f(x) = xexp(— x). Thus we find 
the probability density for a chain segment of length N being formed by breaking a link of strength x 
in the form 

Poo(x\N) = N 2 p(x)P(x) exp[-NP{x)\ . (14) 
Finally, we are in the position to obtain the desired result, the probability density for a segment 
of length N having a weakest link of strength x. This is obtained from the distribution p ao (x\N) as 
follows 

p*(x\N)=J dx -j^ P c (x\x ,N) Poo {x \N) 

= Np(x) f dx exp ( - N[P(x) - P{x )}) Poo (x \N) . (15) 
Jo 

In the limit of large N we find the simple, asymptotically exact result |QJ. 

In the case where p{x) has a finite limit for very weak chains, po — p(0), the strength distribution 
iJTjl has a simple universal scaling form when N 3> 1: 

p*(x\N) ~ ^N 3 plx 2 exp(-Np x) . (16) 
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